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How the vortex lattice of a superconductor becomes disordered: 
a study by scanning tnnneling spectroscopy 

M. Zehetmayer^ 

^ Atominstitut, Vienna University of Technology, 1020 Vienna, Austri^ 

Order-disorder transitions take place in many physical systems, but observing them in detail 
in real materials is difficult. In two- or quasi-two-dimensional systems, the transition has been 
studied by computer simulations and experimentally in electron sheets, dusty plasmas, colloidal 
and other systems. Here I show the different stages of defect formation in the vortex lattice of a 
superconductor while it undergoes an order-disorder transition by presenting real-space images of 
the lattice from scanning tunneling spectroscopy. When the system evolves from the ordered to the 
disordered state, the predominant kind of defect changes from dislocation pairs to single dislocations, 
and finally to defect clusters forming grain boundaries. Correlation functions indicate a hexatic- 
like state preceding the disordered state. The transition in the microscopic vortex distribution is 
mirrored by the well-known spectacular second peak effect observed in the macroscopic current 
density of the superconductor. 


Though studied for many years, the physics of vor¬ 
tex matter of a superconductor— is far from being 
unraveled. There are basically three competing ener¬ 
gies, namely vortex-vortex repulsion, vortex pinning, and 
thermal fluctuations, that determine the state of vortex 
matter—!^. While vortex-vortex repulsion gives rise to an 
ordered vortex lattice, vortex pinning and thermal fluc¬ 
tuations make the lattice disordered. Changing the mag¬ 
netic field, the temperature or the defect structure mod¬ 
ifies the balance of the competing energies, which evokes 
transitions between different vortex matter states. 

The macroscopic critical current of a superconductor 
may have a second maximum (fishtail effect) at a high 
magnetic field, which is the result of an order-disorder 
transition of the vortex lattice^^— . How this transition 
takes place in detail by the formation of lattice defects is 
a fundamental question that has been addressed by the¬ 
oretical worb^ and computer simulation o^°d^ thus far. I 
have addressed this question by experiment, namely by 
pushing the vortex lattice from the ordered to the disor¬ 
dered state in a controlled way in small steps and study¬ 
ing the corresponding real-space lattice images acquired 
by scanning tunneling spectroscopy. 

Order-disorder transitions have been investigated in 
two- or quasi-two-dimensional systems by computer 
simulation a^^d^ and experimentally in electron sheetsi^, 
dusty plasmaai^, colloidali^^— and other system a^^d^ . 
Most results agree with theory^Sr— predicting that the 
lattice becomes disordered in a two-stage process includ¬ 
ing the intermediate hexatic state. The properties of the 
vortex matter are fundamentally different to these sys¬ 
tems, which makes studying its order-disorder transition 
worthwhile. Indeed, the essential point of this work is to 
show in detail the different stages of the order-disorder 
transition of the vortex matter in real-space, driven by 
the successive generation of dislocation pairs, single dislo¬ 
cations, and finally grain boundaries. Additionally, I will 
analyze the different states, including a possible hexatic- 
like state, by evaluating correlation functions and will es¬ 
tablish a direct proportionality between the macroscopic 


critical current and the microscopic disorder of the lat¬ 
tice from the onset to the maximum of the second peak 
region. 

A vortex (flux-line) lattice is formed in a type-II 
superconductor by applying a magnetic field 
Single vortices have a tube-like shape and are roughly 
aligned along the field orientation. They have a normal¬ 
conducting core with a radius in the nanometer range 
and are surrounded by circulating supercurrents gener¬ 
ating a magnetic flux quantum ($0 — 2.07 x 10“^^ Vs). 
Their density is basically proportional to the magnetic 
field via n ~ where fiQ = Air x 10“^ NA“^. 

The circular currents make any two flux-lines repel each 
other, which results in an ordered hexagonal lattice, in 
which no macroscopic currents over macroscopic ranges 
can be accumulated. But pinning of flux-lines by ma¬ 
terial defects opposes that ordering and may result in 
a disordered lattice. Changing, for instance, the exter¬ 
nal magnetic field leads to a macroscopic gradient of the 
vortex density and thus to a macroscopic current when 
pinning is sufficiently strong^. 


RESULTS AND DISCUSSION 
Experiment 

The experiments were carried out on NbSe 2 single crys¬ 
tals (see methods). The critical current density of the as- 
grown (unirradiated) samples had just one maximum at 
a low magnetic field and declined rapidly with increasing 
field. Introducing material defects, capable of pinning 
vortices, by exposing the samples to neutron irradiation 
made the current density increase and, in addition to the 
low-field maximum, a second current peak appear at high 
fields, called the second peak effect (see figure[T]and Sup¬ 
plementary section I). It should be noted that the sam¬ 
ples have a finite thickness of some 100 fim, along which 
vortices can be bent, making the vortex matter system 
more three-dimensional. Double-sided decoration mea- 
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FIG. 1. Comparison of the critical current density and 
the defect density. The critical current density (solid line) 
and the defect density (symbols) of the vortex lattice of the ir¬ 
radiated (Fn = 3 X 10^^ rn“^) NbSe 2 single crystal are almost 
directly proportional from 0.14 to 0.9 T at 4.2 K. 


surements revealed rather stiff vortices in NbSe^^. 

Mapping the density of states of the sample surface, 
scanning tunneling spectroscopy reveals the locations of 
the vortices (see methods). To get the same history- 
dependent vortex matter state as with the macroscopic 
current measurements, it was necessary to first set the 
temperature to 4.2 K and then apply the magnetic field. 
This procedure leads to a small gradient in the macro¬ 
scopic vortex density proportional to the critical cur¬ 
rent, establishing a corresponding driving force on the 
vortices^. 


The vortex lattice and the critical current density 

Figure [5] shows the vortex lattice of the sample of fig¬ 
ure [1] at several magnetic fields from 0.16 T, which is close 
to the minimum of the current, to 0.50 T, which is deep in 
the second peak region. Here each line connects two adja¬ 
cent vortices, hence each intersection marks the position 
of a vortex. As each vortex has six nearest neighbors in 
a hexagonal lattice, a lattice defect is defined as a vortex 
having more or less than six neighbors and is highlighted 
by a large colored symbol in the figures. 

Figure [2] makes it manifest that the disorder, reflected 
by the lattice defect density, increases with field. Figure[T] 
shows that the lattice defect density is almost directly 
proportional to the macroscopic critical current density 
from about 0.14 to 0.9 T. Such a relation is not a-priori 
expected and was indeed not observed at higher fields. 
The microscopic vortex lattice and the macroscopic cur¬ 
rents were determined in the same sample. 


Defects in the vortex lattice and the order-disorder 
transition 

Different kinds of defecti^ appear in the two- 
dimensional vortex lattice images. A single defect, having 
no defect as a direct neighbor, is called a disclination. 
Two adjacent defects from which one has five and the 
other seven neighbors (shown by a red-blue pair in the 
figures) indicate a dislocation. Two connected disloca¬ 
tions form a dislocation pair. Finally, also larger defect 
clusters show up. 

The most established theoretical explanation of how 
a two-dimensional particle system undergoes an order- 
disorder transition was worked out by Kosterlitz, Thou- 
less, Halperin, Nelson, and Young (KTHNY—^— ). In the 
ordered state, characterized by a quasi-long-range trans¬ 
lational and a long-range orientational order, defects ap¬ 
pear only in the form of dislocation pairs. Pushing the 
system towards the disordered state first generates the 
so-called hexatic state by dissolving the pairs into single 
dislocations, making the translational order short-range. 
The final disordered state is dominated by disclinations, 
making also the orientational order short-range. Alterna¬ 
tively, Chui^i proposed that the dislocations form grain 
boundaries instead of dissolving into disclinations. Both 
scenarios were backed by experiments and simulations. 

Previous measurements^ revealed a rather well ordered 
defect-free vortex lattice in an unirradiated clean NbSe 2 
single crystal at fields where the critical current vanished. 
Magnetic decoration images by Kim et al<^ showed that 
defect-free regions may contain more than 10® vortices 
in clean superconductors. In contrast, several dislocation 
pairs and a small number of single dislocations showed 
up in a sample irradiated to a low neutron fluence (Fn = 
10^^ m“^) at fields between the two peaks where the crit¬ 
ical current was also close to zero (Supplementary section 
II). 


The different stages of the order-disorder transition 

Now I come to the details of figures [H and [31 which 
refer to a higher irradiated (Fn = 3 x 10^^ m“^) sample, 
whose current, shown in figure[Tl does not vanish between 
the peaks. At 0.16 T, a field close to the minimum of the 
critical current, several dislocation pairs and a few single 
dislocations appear. To compare lattices with different 
sizes, I normalized the number of defects to a lattice with 
1000 vortices. This gives 13.5 dislocation pairs and 4.5 
single dislocations (figure|31). No other defects are visible. 
Note that all dislocation pairs form loops with alternat¬ 
ing coordination numbers between adjacent defects, thus 
minimizing the Burgers vector— and the corresponding 
defect energy, and destroying the order merely locally. 
It is reasonable that dislocation pairs dominate at low 
disorder, for they can easily be created just by locally 
displacing the involved vortices. In contrast, single dislo¬ 
cations are topological defects, which cannot be created 
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FIG. 2. The order-disorder transition of the vortex lattice in NbSe 2 at 4.2 K. The vortex lattice of NbSe 2 , shown 
at 0.16, 0.18, 0.19, 0.20, 0.25, and 0.50T (cp. with figure [T}, becomes more disordered with increasing field. Each thin line 
indicates a bond between two adjacent vortices, the large (colored) symbols indicate lattice defects, namely vortices with four 
(o), five (o), seven (•), or eight (•) nearest neighbors. The horizontal bars correspond to a length of 1 ^m. The image boundaries 
were not used for the evaluation. 
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FIG. 3. Histograms of the defect cluster size. The 

panels show the number of defect clusters per 1000 vortices 
as a function of the cluster size of all lattices of figure [2] 


by such simple transformational^, but two are obtained 
by splitting a dislocation pair, or they are introduced 
at the sample edges. Likewise, dissolving a dislocation 
leads to two disclinations. The pairs also dominate in 
the lower irradiated sample (Fn = 10^^ m“^) at fields 
between the two current peaks and also form loops with 
alternating coordination number (Supplementary section 
If), but even here single dislocations exist. Some of these 
single dislocations may be remnants of the more disor¬ 
dered state at low fields, where the first current peak ap¬ 
pears. Note that some vortices may be pinned strongly 
enough to maintain their position during the magnetiza¬ 
tion process. 

Increasing the field to 0.18 T leads to merely minor 
modifications. There are still more pairs (10.4) than sin¬ 
gle dislocations (3.2), but, for the first time, also a larger 
defect cluster, consisting of six vortices, emerges. More¬ 
over, some of the pairs no longer exist as loops but form 
lines that are not closed, leading to rather large Burgers 
vectors and defect energies. This may be the first stage 
of a process that dissolves the pairs into single disloca¬ 
tions. It also becomes apparent that the defects are not 
homogeneously distributed in the lattice but rather form 
chains. 


At 0.19 T the changes are more fundamental. In con¬ 
trast to the low-field results, the number of single dislo¬ 
cations (13) exceeds that of the pairs (5) (figureE]). Most 
single dislocations are located close to another one, indi¬ 
cating that they have just been separated. Additionally, 
larger clusters with up to twelve vortices and one disclina- 
tion, located, however, at the image boundary and hence 
presumably being part of a larger cluster, show up. Apart 
from this disclination, all clusters consist of an even num¬ 
ber of vortices. 

With increasing field, the basic trends, observed at 
0.19 T, carry on. The number of single dislocations and 
pairs does not change considerably, while the size of the 
larger clusters grows. The number of clusters decays 
roughly exponentially with their size (figure [3]). Discli¬ 
nations emerge more frequently though do not become 
dominant. 

Going to 0.2 T slightly increases the critical current, 
though the field is still close to the onset of the second 
peak regime. Here, the formation of grain boundaries, 
consisting of defect clusters, becomes manifest. At 0.25 T 
the critical current has increased by more than a factor 
of two and the disorder has grown accordingly. The de¬ 
fect cluster chains become more bulky though defect-free 
grains are still visible. 

Finally, at 0.5 T, which is deep in the second peak 
region, the lattice is quite disordered. The defect den¬ 
sity is about one third, and for the first time, vortices 
with more than seven or less than five nearest neighbors 
emerge. The defect amalgamations are no longer line¬ 
shaped but rather bulky, yet defect-free vortex grains are 
still present. A further increase of the magnetic field to 
1.6 T basically raises the disorder but does not alter the 
nature of the vortex lattice structure. 

Although the finite sample thickness leads to a more 
three-dimensional character of the vortex lines, some sim¬ 
ilarities with KTHNY theory are striking. In particular I 
observed in accordance with KTHNY theory that in the 
most ordered lattice, dislocation pairs dominate, which 
first form closed loops, so that the order is only locally 
destroyed, and later also become elongated. Then, sin¬ 
gle dislocations, destroying the quasi-long-range trans¬ 
lational order, appear more frequently. But finally, in 
contrast to KTHNY theory, which predicts disclocations 
to dissolve into disclinations, though in agreement with 
Chui’s proposal^^, the dislocations become clustered and 
form a grain boundary structure getting more bulky with 
increasing disorder. 


Correlation functions 

Correlation functions offer additional insight into the 
ordering state of our systems (see methods). The trans¬ 
lational correlation function measures the translational 
symmetry of the system and the bond-orientational cor¬ 
relation function the relative orientation of two bonds be¬ 
tween nearest neighbors as a function of distance. Both 
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FIG. 4. The correlation functions of the vortex lat¬ 
tices. The symbols show the (upper) envelope of the ori¬ 
entational (left panel) and of the translational (right panel) 
correlation functions of the unirradiated sample at 1T and of 
the irradiated sample of figure [2] at different magnetic fields. 
The fit functions (solid lines) follow either a power law (red) 
or an exponential (blue) decay behavior. 


functions change from one for a perfectly ordered to zero 
for a completely disordered state. 

Figure 0] shows these functions for our systems. The 
highest order is found in the defect-free lattice of the 
unirradiated sample. Its slightly decaying though quasi- 
long-range translational order and its virtually constant 
orientational order agree with the properties of a Bragg 
glass^SiSS. In the disordered lattices, that is, at and 
above 0.5 T, both functions decay exponentially for not 
too large distances and have similar correlation lengths 
of 1-2 lattice parameters. This would match the vortex 
glass behavior, though the lattice at 0.5 T is not amor¬ 
phous. The flat long-range behavior reflects the fact that 
the mean orientation of different vortex grains does not 
change significantly, presumably an aftermath of the or¬ 
dered state that had been run through during applying 
the field at a constant temperature. Additionally, the fi¬ 
nite sample thickness may also contribute to the unusual 
long-range behavior. At higher helds the translational 
correlation does not change significantly, while the orien¬ 
tational correlation further declines. 

Near the onset of the second peak, for instance when 
the field changes from 0.18 to 0.20 T, the orientational 
order remains almost unchanged and is quasi long-range, 
whereas the translational order decreases appreciably 
and has a much shorter correlation length in both cases. 
This is what one expects of the hexatic state^i. The 
observation of the hexatic state, associated with the ap¬ 
pearance of free dislocations, was also reported for the 
melting process of the vortex lattice in a superconduct¬ 
ing tungsten-based thin film^. The hexatic phase in a 
three-dimensional vortex lattice was discussed in Ref. 1^ . 


Further discussion 

In conclusion, the order-disorder transition of a super¬ 
conductor’s flux-line lattice follows reasonably well the 
theoretical predictions of the grain boundary scenario. 
While in most cases, the systems are cooled down from 
a disordered state, our system reached the final state by 
increasing the magnetic field at a constant temperature 
and thus always first passed the low-field disordered state 
and then the most ordered state at the onset of the second 
peak. But only this procedure makes a comparison of the 
microscopic vortex lattice with the macroscopic critical 
current useful. Moreover, it is a change in temperature 
that evokes the transition in most experiments, whereas 
it is a change in magnetic field in our case. Increasing 
the magnetic field increases the number of vortices and 
thus reduces the distances between them, as it is the case 
when pressure is applied to a conventional system. 

The experiments confirm that the transition is induced 
mainly by pinning and not by thermal effects, since no 
high-held disordered state, at least below 1.7 T, which 
is 85 per cent of the upper critical held, was observed 
in the unirradiated sample, in which vortex pinning is 
weak. Moreover, the disordered vortex lattice conhgura- 
tions did not change signihcantly with time, as was found 
from measuring the same sample area several times in 
immediate succession. This agrees with the common as¬ 
sumption that the order-disorder transition is a solid-to- 
solid transitional^, at least in low-temperature supercon¬ 
ductors. It was shown theoretically^ii that in samples 
with a very small Ginzburg number, such as in NbSe 2 
(10“® - 10“®), thermal effects are not expected to be¬ 
come relevant at helds lower than some 90 per cent of 
the upper critical held, which is far from our transition 
held. Pinning effects extend the region where thermal 
energy matters, but the predicted shift of the melting 
line was found to be small for NbSe 2 even when pinning 
effects are strong^. 

The assumption of a two-stage nature of the order- 
disorder transition in vortex matter was recently con- 
hrmed by two other scanning tunneling spectroscopy 
studie s®®’®"*^ . In the hrst, Guillamon et al<^ investigated 
a tungsten-based superconducting thin him, in which the 
order-disorder transition was generated by a modulation 
of the sample thickness. No lattice defects were found 
at low magnetic helds. In the hrst stage of the transi¬ 
tion, pair and single dislocations showed up, and in the 
second stage, the dislocation density grew strongly and 
free disclinations were created. The translational and ori¬ 
entational correlation functions revealed good agreement 
with KTHNY theory. 

Ganguli et ali^ analyzed a cobalt-doped NbSe 2 sin¬ 
gle crystal having a second peak effect at high magnetic 
helds. They also observed a proliferation of dislocations 
in the hrst stage of the transition and the appearance 
of single disclinations, driving the lattice into an amor¬ 
phous glass, in the second stage. They showed that the 
increasing and the decreasing held branch led to similar 
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vortex structures, but differences were observed in the 
long-range behavior of the correlation functions at fields 
in the transition region, which confirms the assumption 
suggested above that the deviations of the correlation 
functions from KTHNY theory for long distances may 
be explained by the specific history of the magnetization 
process. Cooling the sample at a constant field led to the 
expected higher disorder in the lattices. 


Though all three studies confirm the two-step nature 
of the transition, a detailed comparison reveals some re¬ 
markable differences, which may be caused by different 
material properties and sample thicknesses, leading to 
a more three- IRe f. l34l and this work) or a more two- 
dimensional (Ref. l33h character of the vortex matter. 
First, in the present work, I found the defect structure 
of the rather ordered lattice to be dominated by disloca¬ 
tion pairs, which agrees with the prediction of KTHNY 
theory. In contrast, merely single dislocations showed 
up in Ref. HI at this stage. Second, I found the defects 
to form a grain boundary structure, getting more bulky 
with increasing field. Such a structure was not reported 
by the othergroups, though several vortex lattice im¬ 
ages of Ref. 1^ might indicate an inhomogeneous defect 
distribution. 


Vortex lattice defects concentrated mainly in grain 
boundaries were reported by Fasano et a h^^i^^ in NbSe 2 . 
These results were, however, obtained at very low mag¬ 
netic fields, namely at 3.6 mT, which was not in the sec¬ 
ond peak regime but close to the first critical current 
peak. Employing magnetic decoration, they were able 
to assess large sample areas and thus to detect vortex 
grains including several hundred vortices in both field- 
cooled and zero-field-cooled experiments. Based on the 
assumption that the field-cooled lattices had been frozen 
in close to the upper critical field in the second peak re¬ 
gion, they concluded that the vortex lattice disorder does 
not correlate with the macroscopic current in the second 
peak region, which is in contrast to the proportionality 
found in this work. However, the proportionality between 
current and lattice disorder was only shown for fields be¬ 
tween the onset and the maximum of this region in this 
work, while the frozen lattice of Ref. HI would refer to the 
high-held or high-temperature end of this region, where 
the correlation can be quite different. A grain bound¬ 
ary structure was also found in simulationsi2iiii^“i2S of 
the two-dimensional vortex matter. Other two- or quasi- 
two-dimensional systems such as electron sheets^^, dusty 
plasmaai^i^, colloidaU^^— and other systemsi^d^ were 
found to become disordered by the same or a similar 
two-stage process, though only some of these systems de¬ 
veloped a grain-boundary structure. 


METHODS 

Samples 

The material used for this study is NbSe 2 , which be¬ 
comes superconducting below some 7.2 K. All experi¬ 
ments were carried out at 4.2 K with the magnetic held 
parallel to the uniaxial axis of the single crystal’s hexag¬ 
onal unit cell, for which the upper critical held is some 
2 T. Recently, NbSe 2 was shown to be a two-band super¬ 
conductor similar to MgB 2 ^. Our samples had a size of 
1-3 mm in lateral direction and 100 - 200/im parallel to 
the uniaxial axis. 


Neutron irradiation 

Neutron irradiation of the samples took place in a 
TRIGA MARK H research hssion reactor—. This proce¬ 
dure introduces crystal defects that are efficient for vor¬ 
tex pinning in many superconductors. For instance, in 
MgB 2 ^ and YBa 2 Cu 307 _,?^, transmission electron mi¬ 
croscopy revealed spherical defects with a diameter close 
to the superconducting coherence length. Neutron ir¬ 
radiation tremendously strengthens vortex pinning also 
in NbSe 2 , as shown by its effects on the critical current 
(Supplementary section I). One sample was irradiated to 
a fast neutron {E > 0.1 MeV) huence of 1 x 10^^ m“^ 
and the sample to which most results of this article refer 
to 3 X 10^^ m~^. 


Critical current density 

The critical current density was obtained from SQUID 
magnetometry by measuring the magnetic moment as a 
function of magnetic held at 4.2 K. Assuming the current 
to be constant in the sample, one can derive the critical 
current from the hysteresis of the magnetic moment in 
increasing and decreasing held^. The unirradiated sam¬ 
ples showed a hnite though small critical current only 
at low helds. A low neutron huence enhanced the cur¬ 
rent signihcantly at low helds and additionally produced 
a second peak close to the upper critical held. With in¬ 
creasing huence the low-held current became larger and 
the second peak more prominent. 


Scanning tunneling spectroscopy 

Scanning tunneling spectroscopy reveals the electronic 
density of states at the sample surface. In a superconduc¬ 
tor the procedure produces the well-known energy gap of 
the quasi-particles near the Fermi level. This gap, how¬ 
ever, does not exist in the vortex core, and hence the 
tunneling signal changes, thus identifying the vortex core 
positions, which makes scanning tunneling spectroscopy 
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the sole known method for measuring the vortex lattice 
in real-space at arbitrary magnetic fields^i^i^^— . 

For the measurements, a sample was cleaved in air and 
then put into the sample tube, which was first evacuated 
and then flooded with pure helium gas. The measure¬ 
ments were carried out with a Ptir tip. Having placed 
the sample in a helium cryostat, I cooled it to 4.2 K and 
then applied a magnetic field. Then, after waiting for 
several hours, the surface topography was mapped by a 
constant-current mode measurement and, at the same 
time, the density of states at an energy level slightly 
above the Fermi level but below the energy gap width 
of some 2 meV by a lock-in technique (Supplementary 
section III). The densities of state maps usually showed 
a clear contrast for the vortex core positions (Supple¬ 
mentary section III). Applying simple filters such as a 
Gaussian filter to the raw data reduced the noise, so that 
the vortex core positions were easily acquired by seeking 
all local maxima. The nearest neighbors of each vortex 
were determined by a Delaunay triangulation. 

Correlation functions 

The translational correlation functional is defined as 
GT(|r-f^|) = (e*^“e-'^“') (1) 


where G is a reciprocal lattice vector, and u is the dis¬ 
placement of a vortex at r from its lattice position i?, 
that is, u = f — R and = 1. The brackets indi¬ 
cate an averaging over all vortex pairs. The reciprocal 
lattice vector was determined by seeking the vector with 
the largest value near the boundary of the first Brillouin 
zone after a Fourier transformation of the lattice. The 
distances were discretized into interval lengths of 0.1 oq. 

For the bond orientational correlation functionli, one 
has to calculate 


Ge(|F-r-|) = (66(r)6^(r')) (2) 


with 


1 ” 

6e(r) = -V (3) 

n 


The sum runs over all nearest neighbors j of the vortex 
at r, and r)) is the angle of the bond connecting the 
two vortices with respect to the x-axis. 
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